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Cosmology is making impressive progress and it is producing stringent bounds on the sum of the 
neutrino masses E, a parameter of great importance for the current laboratory experiments. In 
this letter, we exploit the potential relevance of the analysis of Palanque-Delabrouille et al. [JCAP 
1502, 045 (2015)] to the neutrinoless double beta decay (Ov/Iff) search. This analysis indicates small 
values for the lightest neutrino mass, since the authors find E < 84meV at lcrC.L., and provides a 
la preference for the normal hierarchy. The allowed values for the Majorana effective mass, probed 
by 0 vfdf), turn out to be < 75meV at 3 ctC. L. and lower down to less than 20meV at IctC.L.. 

If this indication is confirmed, the impact on the 0 is/3/3 experiments will be tremendous since the 
possibility of detecting a signal will be out of the reach of the next generation of experiments. 
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I. INTRODUCTION 


Neutrinoless double beta decay (Oi'/J/J) [1] probes lep¬ 
ton number conservation and allows the investigation of 
the nature of the neutrino mass eigenstates. Assuming 
that “ordinary” light neutrinos mediate the 0^/3 /3 transi¬ 
tion, the key parameter that regulates the decay rate is 
the Majorana Effective Mass, namely 




e iai |f/ei|mi + e M2 |C/ e 2 2 |m 2 + |t/ 2 3 |m 3 


( 1 ) 


It represents the absolute value of the ee entry of 
the neutrino mass matrix. Here, m,; are the masses 
of the individual neutrinos Vi, ai i2 are the Majorana 
phases and U e i are the elements of the mixing matrix 
that defines the composition of the electron neutrino: 
\v e ) = From Eq. (1), it is evident that 

the 0 v/3/3 amplitude strongly depends upon the Majorana 
phase variation and the absolute neutrino mass scale. 

The mixing parameters and the squared mass dif¬ 
ferences are effectively measured in oscillation experi¬ 
ments [2], Since a\ and a 2 are very challenging to test [3], 
we focus on the implications coming from measurements 
of neutrino masses which, at present, are best probed by 
cosmological surveys. 


II. RESULTS FROM COSMOLOGICAL DATA 

The three light neutrino scenario is consistent with all 
known facts in particle physics including the new mea¬ 
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surements by Planck [4]. In this assumption, the quantity 
E is just the sum of the masses of the neutrinos: 

E = mi + m 2 + m 3 . (2) 

The present information on three-flavor neutrino oscil¬ 
lations is compatible with two different neutrino mass 
spectra: normal hierarchy ( J\fTL ) and inverted hierarchy 
(TH). In the former case one gets: 

! toi = m 

to 2 = \/ m 2 + 5m 2 (3) 

to 3 = V m 2 + Am 2 + 5m 2 /2 


while, in the latter: 


{ mi = y/m 2 + Am 2 — 5m 2 / 2 
to 2 = \Jm 2 + Am 2 + 5m 2 / 2 (4) 

m 3 = m. 

Here, m is the lightest neutrino mass, while 5m 2 and 
A to 2 are the mass splittings measured by oscillations, 
defined according to Ref. [2]. 

A recent analysis [5], based exclusively on oscillation 
data, finds a lcr preference for the 1TL. This conclu¬ 
sion relies on partially contradictory information. In fact, 
the analysis of atmospheric neutrinos performed by the 
Super-Kamiokande Collaboration [6] shows a “tendency” 
(i.e. a lcr preference) in favor of the NIL. 

The determination of the mass hierarchy is a difficult 
challenge. In view of its great importance and of the 
following discussion, we believe that future global anal¬ 
yses could benefit from the inclusion of the information 
obtained from cosmology. 

In this work, we focus on the tightest experimental lim¬ 
its currently available for E. These are usually obtained 
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by combining cosmic microwave background (CMB) data 
with the Lyman-a forest ones. Since they probe different 
length scales, the combination allows a more effective in¬ 
vestigation of the neutrino induced suppression in terms 
of matter power spectrum, both in scale and redshift. 

Among the different analyses performed, it is worth 
stressing the following 2a C. L. upper limits: 0.17 eV ob¬ 
tained in Ref. [7] by combining CMB data of the Wilkin¬ 
son Microwave Anisotropy Probe, galaxy clustering and 
the Lyman-a forest of the Sloan Digital Sky Survey 
(SDSS); 0.18eV of Ref. [8] using Planck and WiggleZ 
galaxy clustering data; 0.14 eV obtained in Ref. [9] by 
combining Lyman-a SDSS data with Planck; 0.17 eV ob¬ 
tained in Ref. [4] by using Planck temperature and po¬ 
larization measurements including a prior on the Hubble 
parameter, Supernovae and Baryonic Acoustic Oscilla¬ 
tions (BAOs). 

More recently, by using a new sample of quasar spec¬ 
tra from SDSS-III and Baryon Oscillation Spectroscopic 
Survey searches and a novel theoretical framework which 
incorporates neutrino non-linearities self consistently, 
Palanque-Delabrouille et al. [10] have obtained a new 
tight limit on E. This constraint was derived both in fre- 
quentist and bayesian statistics by combining the Planck 
2013 results [11] with the one-dimensional flux power 
spectrum measurement of the Lyman-a forest of Ref. 

I 12 ]- 

The one-dimensional flux power spectrum is the key 
observable used to derive the tight upper bound on E. 
The final results do not depend upon cosmological and 
astrophysical parameters since the constraints obtained 
on the total neutrino mass E derived in Ref. [10] include 
both statistical and systematic uncertainties. The error 
budget is dominated by systematic uncertainties, that 
are treated as nuisance parameters affecting the observed 
flux power spectrum. 

A detailed modeling and an accurate description of the 
sources of uncertainties that affect this measurement is 
quantitatively presented in Ref. [12] and we refer to this 
paper for a more exhaustive discussion (in particular, see 
Tables 7 and 8 in there). In this paper, we just highlight 
the crucial aspects taken under consideration, such as 

i ) the continuum fitting that allows to remove the 
long wavelength dependence induced by the quasar, 

ii) the metal contamination, 

Hi) the thermal history of the intergalactic medium, 

iv ) the instrumental resolution and noise, 

v ) the impact of galactic feedback (Active Galactic 
Nuclei and Supernovae) and fluctuations of the ul¬ 
traviolet background, 

vi ) the impact of strong absorbers, 

vii) the numerical modeling of small scales in hydrody- 
namical simulations. 


As a final remark, Ref. [12] compares two different meth¬ 
ods of power spectrum extraction and contrasts with the 
previous results by Ref. [13], finding very good agree¬ 
ment. For a comprehensive assessment of astrophysical 
effects on the flux power, we refer also to other analy¬ 
ses in Refs. [10, 14-18], while for uncertainties related 
to the method used in obtaining the limits, we refer to 
Refs. [10, 16]. 

For the purpose of the present paper, it is important 
to stress that these effects have been modeled by care¬ 
fully quantifying their impact in terms of flux power from 
state-of-the-art hydrodynamical simulations. The ampli¬ 
tude and shape of such effect is then parameterized (tem¬ 
plate fitting) and marginalized over in the likelihood cal¬ 
culation. As a result, all the quoted relevant numbers 
already include the effect of the systematic uncertainties 
and the dependence of the results upon the other cosmo¬ 
logical parameters. The relatively wide range of scales 
and redshifts spanned by the observed one-dimensional 
flux power spectrum allows to effectively break the de¬ 
generacies among astrophysical and cosmological param¬ 
eters, providing tight constraints in terms of neutrino 
masses. 

We would like to comment on the widespread attitude 
to dismiss cosmological measurements, invoking the pos¬ 
sible existence of other sources of systematics, not yet 
identified, that could affect to some extent the conclu¬ 
sion. This is in principle possible, as it is for any kind 
of measurement. In some sense, it can be considered 
a super-conservative approach. However, we believe that 
the best thing that can be done is identifying the possible 
relevant systematics and to quantify their impact, just as 
it was done in Ref. [12]. This is especially important to 
continuously improve the reliability of the present cos¬ 
mological measurements, also in view of their impact on 
the interpretation of Ojz/ 3/3 in terms of massive neutrinos. 

From the frequentist interpretation (which is in excel¬ 
lent agreement with the bayesian results), the authors of 
Ref. [10] compute a probability for E that can be sum¬ 
marized in a very a good approximation by: 


Ay 2 (E) 


(E — 22 meV) 2 
(62 meV) 2 


(5) 


Starting from the likelihood function £ oc exp—(Ay 2 /2) 
with Ay 2 as derived from Fig. 7 of Ref. [10], one can 
obtain the following limits: 

E < 84 meV (lcrC.L.) 

E < 146 meV (2trC.L.) (6) 

E < 208 meV (3crC.L.) 


which are very close to those predicted by the Gaussian 
Ay 2 of Eq. (5). 

It is worth noting that, even if this measurement is 
compatible with zero at less than lcr, the best fit value is 
different from zero, as expected from the oscillation data 
and as evidenced by Eq. (5). 
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FIG. 1. (Left) Allowed regions for mpp as a function of S with constraints given by the oscillation parameters. The darker 
regions show the spread induced by Majorana phase variations, while the light shaded areas correspond to the 3cr regions due to 
error propagation of the uncertainties on the oscillation parameters. (Right) Constraints from cosmological surveys are added 
to those from oscillations. Different C.L. contours are shown for both hierarchies. Notice that the lcr region for the XT-i case 
is not present, being the scenario disfavored at this confidence level. The dashed band signifies the 95% C.L. excluded region 
coming from Ref. [10]. 


Furthermore, the (atmospheric) mass splitting A = 
V Am 2 cs 49 meV [2] becomes the dominant term of Eqs. 
3 and 4 in the limit m —i 0. Under this assumption, 
in the case of AfH (XH) E reduces approximately to A 
(2A). This explains why this result favors the A fH mass 
spectrum, as pointed out in Ref. [10] and as advocated 
in older theoretical works [19]. 

It is the first time that some data probing the absolute 
neutrino mass scale indicate a preference for one spe¬ 
cific mass hierarchy. Nonetheless, these results from cos¬ 
mology have to be taken with due caution. In fact, the 
scientific literature contains several authoritative claims 
for a non-zero value for E, see e.g. Refs. [20-23]. In 
particular, it has been very recently suggested [22, 23] 
that a total non zero neutrino mass around 0.3 eV could 
alleviate some tensions present between cluster number 
counts (selected both in X-ray and by Sunyaev-Zeldovich 
effect) and weak lensing data. In some cases, a ster¬ 
ile neutrino particle with mass in a similar range is also 
advocated [11, 24]. However, evidence for non-zero neu¬ 
trino masses either in the active or sterile sectors seems 
to be claimed in order to fix the significant tensions be¬ 
tween different data sets (CMB and BAOs on one side 
and weak lensing, cluster number counts and high val¬ 
ues of the Hubble parameter on the other). In fact, ex¬ 
tending the model by including massive neutrinos does 
not improve much the agreement among different observ¬ 
ables, since CMB and BAOs data do not support this 
extension [25]. This could suggest that systematic errors 
are not fully under control. More precise measurements 
from cosmological surveys are expected in the near fu¬ 
ture (among the others, DESI 1 and the Euclid satellite 2 ) 


and they will probably allow more accurate statements 
on neutrino masses. 


III. IMPLICATION FOR THE Ou/3/3 SEARCH 

The close connection between the results on neutrino 
mass obtained in the laboratory and those coming from 
cosmology was outlined long ago [26]. In order to dis¬ 
cuss the mass parameter mpp which is relevant for the 
0p/3p, it is useful to adopt the representation introduced 
in Ref. [27], where mpp is expressed as a function of E. 
The result is shown in the left panel of Fig. 1. The pro¬ 
cedure to obtain this plot is quite straightforward: 

- the starting points are the expressions for the max¬ 
imum and the minimum values of mpp as a func¬ 
tion of the lightest neutrino mass to, the oscillation 
parameters and the unknown Majorana phases. 
These last are left free to vary between their ex¬ 
tremal values, thus providing an interval of values 
for mpp for a given set of to and oscillation pa¬ 
rameters (see e. g. Refs. [3, 28] for a more detailed 
explanation). 

- for both the mass hierarchies, one has to solve the 
quartic equation that gives m as a function of E, 
Sm 2 and Am 2 . This solution is replaced in the 
previous expressions of mpp, which in turn become 
a function of E and of the oscillations parameters. 

- the residual uncertainties on the oscillation param¬ 
eters must be propagated on the maximum and 
minimum values of mpp (the procedure is described 
e.g. in Ref. [3]). This operation slightly widens the 
allowed regions for mpp. 


1 http://desi.lbl.gov 

2 http://www.euclid-ec.org 
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Mass spectrum mpp max 

[meV] (C. L. 

on E) 


la 

2cr 

3 a 

a m 

16 

41 

64 

XU 

- 

57 

75 


TABLE I. Maximum values for mpp once the new constraints 
on E from Ref. [10] are added. The lcr C. L. maximum value 
on mpp for the XXL is not reported because the scenario is 
excluded in this case. 


constant [29-31], which would imply even longer lifetimes 
and conversely would worsen the sensitivity to mpp sig¬ 
nificantly [3]. 

The impact of the new constraints on E appears even 
more evident by plotting mpp as a function of the mass of 
the lightest neutrino [28]. In this case, Eq. ( 8 ) becomes: 

(r%W ) 2 rn 2 
(na[mpp(m)]) 2 m(E „) 2 


- the minimum value of E, namely E(m = 
0, <5m 2 , Am 2 ), is different from zero for both the 
mass hierarchies and its value is subject to the 
residual uncertainties on 5m 2 and Am 2 . These can 
be easily propagated as follows: 



\d5m 2 


a {5m 2 ) 


f d E 


V d Am 2 


a(Am 2 ) ) . (7) 


This explains the widening of the shaded areas in 
the leftmost part of the allowed regions in the left 
panel of Fig. 1. 

The new cosmological constraints on E from Ref. [10] 
can be now included at any desired confidence level by 
considering the following inequality: 

(y-m^(S )) 2 (E-E(O )) 2 

(M[m w (S )]) 2 + (E n -E(0 )) 2 ^ 

where m/ 3 ^(E) is the Majorana Effective Mass as a func¬ 
tion of E and a[mpp( E)] is the lcr associated error, com¬ 
puted as discussed in Ref. [3]. E„ is the limit on E de¬ 
rived from Eq. (5) for the C. L. n = 1, 2, 3,... By solving 
Eq. ( 8 ) for y , it is thus possible to get the allowed contour 
for mpp considering both the constraints from oscillations 
and from cosmology. In particular, the Majorana phases 
are taken into account by computing y along the two ex¬ 
tremes of mpf 3 (E), namely and m^ n (E), and 

then connecting the two contours. The resulting plot is 
shown in the right panel of Fig. 1. 

The most evident feature of Fig. 1 is the clear differ¬ 
ence in terms of expectations for both mpp and E in 
the two hierarchy cases. The relevant oscillation param¬ 
eters (mixing angles and mass splittings) are well known 
and they induce only minor uncertainties on the expected 
value of mpp. These uncertainties widen the allowed con¬ 
tours in the upper, lower and left sides of the picture. 
The boundaries in the rightmost regions are due to the 
new information from cosmology and are cut at various 
confidence levels. It is notable that at lcr, due to the 
exclusion of the XXL, the set of plausible values of mpp 
is highly restricted. Table I summarizes the new allowed 
values for mpp. 

The next generation of Oir/3/3 experiments is expected 
to probe the upper values of the predicted XXL region and 
reach a sensitivity for mpp of about 70meV [3]. This as¬ 
sumes the absence of the quenching of the axial coupling 


Here, mpp(m) is mpp expressed as a function of m, 
<j[mpp(m)\ is the lcr associated error, computed as dis¬ 
cussed in Ref. [3], and m(E„) is the value of m calcu¬ 
lated for a given E„. In Fig. 2 the expected sensitiv¬ 
ities for two examples of next generation 0 is/3/3 experi¬ 
ments (CUORE [32] and GERDA-II [33]) are presented. 
The dashed contours indicate the 3cr regions allowed by 
oscillations. The shaded areas are the new allowed re¬ 
gions once the cosmological constraints are added. The 
plot globally shows that the next generation of experi¬ 
ments will have small possibilities of detecting a signal of 
due to light Majorana neutrino exchange. There¬ 
fore, if the new results from cosmology are confirmed or 
improved, ton or even multi-ton scale detectors will be 
needed [3]. 

On the other hand, a 0signal in the near future 
could either disprove some assumptions of the present 
cosmological models, or suggest that a different mech¬ 
anism other than the light neutrino exchange mediates 
the transition. New experiments are interested in test¬ 
ing the latter possibility by probing scenarios beyond the 
Standard Model [34-36]. 



FIG. 2. Constraints from cosmological surveys are added to 
those from oscillations in the representation mpp as a func¬ 
tion of the lightest neutrino mass. The dotted contours repre¬ 
sent the 3 <j regions allowed considering oscillations only. The 
shaded areas show the effect of the inclusion of cosmologi¬ 
cal constraints at different C. L. . The horizontal bands cor¬ 
respond to the expected sensitivity for future experiments, 
according to Ref. [3]. 
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IV. CONCLUSION 

A cautious attitude in dealing with the results from 
cosmological surveys is highly advisable. However, the 
newest analysis [16] confirms and strengthens the cosmo¬ 
logical indications, and it is likely that we will have soon 
other substantial progress. Moreover, the present theo¬ 
retical understanding of neutrino masses does not contra¬ 
dict these cosmological indications. These considerations 
emphasise the importance of exploring the issue of mass 
hierarchy in laboratory experiments and with cosmolog¬ 
ical surveys. 

From the point of view of Oiz/3/3, these results show that 
ton or multi-ton scale detectors will be needed in order 
to probe the range of rnpg now allowed by cosmology. 
Nevertheless, if next generation experiments see a signal, 
it will likely be a 0 vj3f3 signal of new physics different 
from the light Majorana neutrino exchange. 
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Note added: After this work was completed, new com¬ 
pletely independent analysis with results fully consistent 
with the one discussed in Ref. [10] were presented. In¬ 
deed, in Refs. [37] and [38], the limits E < 0.11 and 
E < 0.113 eV at 95%C.L., respectevely, were derived by 
combining data from BAOs, CMB and galaxy clustering. 
This adds confidence to our hypotheses and inferences. 
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